Testing two nuclear physics approximations used in the standard 
leaky box model for the spallogenic production of LiBeB 



J.P. Kneller^ 
kneller@pacific.mps . ohio-state . edu 
g : and 

o 

(N ; J.R. Phillips^ 

X5 \ 

! phillips@pacific.mps.ohio-state.edu 

^ : and 

> ■ T.P. Walkeri'2 

^ . 

Q , twalker@pacific.mps.ohio-state.edu 

(N 
O 
cn 
o 



ABSTRACT 



p i' 

I ' The spallative production rates of Lithium, Beryllium and Boron (LiBeB) are 

^ ■ a necessary component in any calculation of the evolution of these nuclei in the 

^ ■ Galaxy. Previous calculations of these rates relied on two assumptions relating to 
the nuclear physics aspects: the straight-ahead approximation that describes the 

yj ■ distribution of fragment energies and the assumption that the major contributor 

j_j ■ to the production rate arises from single-step reactions between primary cosmic 



ray projectiles and interstellar medium targets. We examine both assumptions 
by using a semi-empirical description for the spall's energy distribution and by 
including the reactions that proceed via intermediary fragments. After relaxing 
the straight-ahead approximation we find the changes in the production rates and 
emerging fluxes are small and do not warrant rejection of this approximation. 
In contrast we discover that two-step reactions can alter the production rate 
considerably leading to noticeable increases in the efficiency of producing the 
LiBeB nuclei. Motivated by this result we introduce a cascade technique to 
compute the production rates exactly and find that the results differ only slightly 
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from those of our two-step calculations. We thus conclude that terminating the 
reaction network at the two-step order is sufficiently accurate for current studies 
of spallation. 

Subject headings: abundances — cosmic rays — nuclear reactions 

1. Introduction 

The production of the isotopes of Lithium, Beryllium and Boron (LiBeB) from the 
spallation of ^^C, ^^N and ^^O in the Inter-Stellar Medium (ISM) by protons and alpha par- 
ticles in the Cosmic Rays (CRs) outlined in Meneguzzi, Audouze & Reeves (1971) filled the 
last major gap in the understanding of the origin the elements. This method of manufac- 
ture (dubbed Galactic Cosmic Ray Nucleosynthesis (GCRN) or Cosmic Ray Spallation) is 
thought to be the sole source of ^Li, ^Be and ^'^B and thus observations of their abundances 
give clues to the contribution of spallogenic ^Li and ^^B to the total abundance of these 
nuclei. The fragihty and paucity of these elements renders the derivation of their abun- 
dances from stellar spectra challenging but despite this difficulty data now exists over a wide 
range of metallicity (Molaro et al. 1997; Duncan, Lambert & Lemkc 1992; Duncan et al. 
1997; Garcia Lopez et al. 1998; Hobbs, Thorburn & Rcbull 1999; Nisscn et al. 1999; Smith, 
Lambert & Nissen 1993). After comparing the observations with the predictions of theory 
many conclusions have been drawn: for example the discrepancy between the observed solar 
^^B/^'^B isotopic ratio of 4.05 (Anders & Grevesse 1989) with the spallogenic prediction of 
~ 2.5 (Reeves 1974) indicates another source of ^^B must exist. 

The strength of this, and indeed any, conclusion ultimately relies on the validity of the 
underlying approximations that permeate the GCRN calculation of which there are many. 
Within the sources of error there are two that emerge from the nuclear physics aspects of 
the problem. The first is the straight- ahead approximation which assumes that the fragment 
emerges from a reaction with a velocity equal to that of the projectile, while the second is 
that the only significant source of LiBeB is the direct, one-step, production from the spalling 
of the CNO and that production via intermediaries is insignificant. The use of both of these 
approximations is widespread but remains largely untested. Relaxing either considerably 
increases the complexity of the calculation yet this is exactly our intention in this paper. In 
section §2 we present the basics of cosmic ray calculations in the leaky box model, show how 
the reaction expansion is invoked and how the straight-ahead approximation enters. Then in 
section §3 we relax the straight-ahead approximation and demonstrate the implications for 
the production rates and fiux of the secondary nuclei. In section §4 we include the two-step 
processes in the reaction rates and then, motivated by these results, we proceed to section 
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§5 where we use a cascade technique to calculate the production rates exactly. 



2. CR Basics and the Straight-Ahead Approximation 



The theoretical underpinning for the study of cosmic rays is a propagation model that 
relates a source spectrum to the flux observed at a later time and place. Of the theoretical 
models used the most popular for the calculation of the rates of production of ^Li, ^Li, ^Be, 
^°B and ^^B is the Leaky-Box Model (see e.g. Cesarsky (1980)) or its variants. The model 
is simpler and more transparent than the rival diffusion model but it suffers from number of 
deficiencies, both theoretical and when used to fit some cosmic ray data. Despite this failing 
the use of the model is widespread because it has empirically proven to be useful in fitting 
compositional cosmic ray data which is the focus of this paper. 

In the leaky box model the flux of a cosmic ray species is given by (Fields, Olive and 
Schramm 1994; Ramaty et al. 1997) 



where Q represents the source of particles, the mass density of the interstellar medium 
(ISM) is represented by p and the integrating factor S has a physical interpretation as the 
probability that a particle produced with energy E will survive to form part of the ISM. 
This survival probability is simply 



where w is the energy per nucleon lost per unit distance traversed and per unit mass density 
of the medium and A(£') is the path length or grammage. The quantity (pw p is frequently 
called the current, the rate of flow of particles from high to low energy, so that in the limit 
when S{E) — * 1 all the particles at a given energy will eventually accumulate at E = i.e. 
they form part of the ISM. This limit is achieved when £^ — > so the rate of change of the 
abundance y is therefore 



where uh is the number density of hydrogen in the ISM. 

So far we have not distinguished between the different nuclei under consideration since 
equations (1), (2) and (3) are valid for each regardless of how the different species become 
cosmic rays. There are differences of course, the flux of CNO is almost entirely that accel- 
erated from the ISM whereas the flux of LiBeB is almost entirely produced 'in-flight' so to 
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(2) 




(3) 
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speak. We can thus make a distinction between two types of sources Q: a primary spectrum 
from whatever mechanism is responsible for acceleration of particles from the ISM, and 
a secondary production term Q^*-* that is due to particle reactions between primary CRs 
and the constituents of the ISM. Cosmic rays are classified by whichever is dominant. The 
production term Q*-*-* is simply the sum of all the reactions that can make the nucleus we 
are interested in, so we have 

Qf{Ef) = Qf{EF)-rQf{EFl (4) 

/"OO 

Q^fiiEp) = dEpMEp)nTap^{Ep)Pp{EF\Ep) (5) 

P,T ''^ 

where we have used the superscripts/subscript P, T and F to denote the quantities associated 
with the projectile, target or fragment of the reaction P + T — > F + X, nr is the number 
density of targets, trp^ is the cross-section for the reaction and Pf{Ef\Ep) is the probability 
distribution for the fragment energy Ep at a given projectile energy Ep. At this point we 
expanded the projectile flux (f)p as (f)p — X]i=o^P' reasoning will become clear in a 
moment, and insert this expression into equation (5) so that the secondary source for the 
fragment Qp^ becomes 

POO 

Q^;\Ep) = dEp(t>f{Ep)nTa^/{Ep)P{Ep\Ep) = Y,Q?, (6) 

i=0 P,T i=l 

/"OO 

QfiEp) = dEp4-'\Ep)nTa^/{Ep)P{EF\Ep), (7) 

P,T 

which introduces the expansion over Q^p . Note that this expansion begins aXi — 1 not i = 
because Q^P is already in use as the primary source spectrum for nucleus F. We can now 
bundle together the primary and secondary sources we introduced in equation (4) into one 
expression: Qp = Xli=o Qf''- Now we can also expand the fragment's flux (pp in exactly the 
same fashion as the projectile flux (pp so that equation (1), now reads 

MEP) = Y.^'i^{Ep) = Y: HdE'pQfiE'p) (8) 

If we equate terms in i we achieve our final result 



= ^TTFT^ / dE'pQ^FiE'p).^^^ (9) 



The meaning of the expansion now becomes clear, we have expanded the fragment fiux (pp 
over the number of reactions that lead from a given projectile P to the desired fragment F 
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because Q^p , defined in equation (7), depends upon the fiuxes of projectiles (f)p~^^. Thus (f)p^ 
is the flux that is accelerated from the ISM, (p^P denotes the flux produced by interactions 
involving (j)p\ (j)p^ is that which arises from reactions involving (f)p^ and so on. 

Inserting the same expansion of into equation (3) leads to a similar expansion of 
the rate of change of abundance dyp/dt as a sum of terms dy^p /dt where each contribution 
is from the source term Q^p 

riE-,Qf(E',)SHE'.). (10) 

We return our attention to equation (7) and the first term in the reaction source expan- 
sion Q^P . This is the term that corresponds to the formation of F via a single interaction of 
a primary CR and the ISM, i.e. the one-step source. We can break this quantity apart into 
the sum from each projectile/target pair and likewise expand (f)p^ and dyp^ /dt in the same 
fashion with each term representing the contribution from the source term Qp . Lastly, the 
primary projectile flux ^p-* may be rewritten as ^p"* = apipp so that the integral over the 
projectile's energy per nucleon Ep yields ap, the total number of projectiles per unit time 
and per unit area. The parameter ap acts as the flux strength and we can therefore remove 
the trivial dependence upon the chemistry of the CRs. The rate of change of the abundance 
in the ISM and the flux of nucleus F are thus 

, (1) , P,T 

^Vf st^ ^Vf sr^ T-iP.T /-,-,\ 



where 



/•oo f poo ~\ 

R^/ ^ J dE'plj dEpMEp)4'^iEp)Pp{E'p\Ep)i SpiE'p) (13) 
^r{Ep) = HdE'p I r dEpMEp)<7p'^{Ep)Pp{E'p\Ep)\ J U) 

J Bp I Jo J Wf[Ef) bF[EF) 



Equations (13) and (14) serve to deflne the normalized, one-step reaction rate Rp^ and the 



normalized, one-step flux "ip 

With the basics complete we can flU in some of the mundane details of the calculation. 
The set of nuclei we will use for both the ISM and CRs are p/^H and a/^He plus all the 
stable nuclei from ^Li through to ^^O. In addition we must add ^Be and ^^C to our set of 
nuclei found in the CRs because their decay mode is inner orbital electron capture and so 



-6- 



are stable in the ionized state. We will sometimes denote by Z the subset of nuclei with 
mass greater than ^He. The stopping power wf{E) is tabulated up to 12 MeV/nucleon for 
each fragment by Northcliffe & Shilling (1970). Above this energy we use the stopping power 
for protons as tabulated by Janni (1982) and rescale in exactly the same manner as Fields, 
Olive and Schramm (1994). The important cross-sections needed for our calculation were 
tabulated by Read & Viola (1984) but there is no cross-sectional data involving ^Be and ^^C 
as a reactant so we have appealed to the fact that each is the mirror of a nucleus for which 
data does exist and therefore assume that the cross-section for reactions involving ^Be or 
^^C is the same as for it's isobaric brother. For the primary CR spectrum we adopt the form 
presented by Gloeckler & Jokipii (1967) and used by Walker, Mathews & Viola (1985) and 
Steigman & Walker (1992) 

c^^p\Ep) = ap— 1:^^^ (MeV/nucleon)-^ cm-^s-^ (15) 

We have set the parameter Eq to be the miclcon mass ttin though we note Garcia-Munoz et 
al. (1987) preferred a value of 400 McV/nuclcon. 

The last quantity we need to specify is the fragment energy distribution Pp{Ep\Ep). 
This quantity has been historically approximated with a 5-function form known as the 
straight-ahead approximation. The approximation is divided into two, the selection being 
dependent upon the identity of the projectile P, 

5{Ep) PG{p,a}1 

5{Ep -Ep) Ve{L] ^ ' 

Note that in the case of P e {p, a} there is no flux of the fragment F. 

We compile the one-step production rates as calculated by equation (13) using 

two, constant values of A = 5 g/cm^ and A = 10 g/cm^. In table (1) only those where the 
reactants are p/^H and ^^O. An initial examination of the results shows that loss of the 
fragments when P e {Z} reduces the production rate by a factor < 10 relative to the rate 
when the projectile and target are interchanged. While this may appear to render the inverse 
reactions unimportant it must be remembered that, at the present time, the abundances of 
many of the Z nuclei in the CR are enhanced relative to their ISM values by approximately 
equivalent factors. Finally, we also see from the table that increasing the grammage leads to 
a larger fraction of LiBeB captured from the inverse reactions and so increasing the efficiency 
of production. 



Pp(yEp\Ep) 



- 7- 



3. Relctxing The Straight Ahead Approximation 

In section §2 we used the straight-ahead approximation to calculate Rp and \E'y 
and we now wish to test its validity. An examination of the approximation was made by 
Fields, Olive and Schramm (1994) who replaced it with another ^-function but with different 
relation between the fragment and primary energies while Tsao et al. (1995) conducted a 
test for heavy nuclei that included a spread in fragment energy. It is the direction of Tsao 
et al. (1995) that we shall explore. 

The momentum distributions of the isotopes produced by fragmentation of various pro- 
jectile and target nuclei are discussed by Morrissey (1989) and Hufner (1985). The distribu- 
tions from experiment are found to be normally distributed in the rest frame of the spalled 
nucleus Z with narrow dispersions and a small mean momentum in the direction parallel to 

the incident projectile. In particular the momentum distributions of the fragments from -"^^C 
and ^^O projectiles upon various targets ranging in mass from Be to Pb were measured by 
Greincr eA, al. (1975) who found that their results had no significant correlation with target 
mass or beam energy. Goldhaber (1974), and more recently Bauer (1989), explain these 
results in terms of a nuclear model with minimal correlations between the nucleon momenta. 
They predict the dependence of the momentum per nucleon distribution Pf{P^p\Ep) of a 
fragment with mass Ap to be 

PF{^f\EF) = --i^exp ( - ^P^^' 7 P^^' 1 (17) 



(2 7re|)3/2-^l 2e 



(z) ~ 
where Y* p is the momentum per nucleon of the fragment in the rest-frame of Z, Pp is the 

mean and e\ the variance of the distribution. To obtain Pf{Ep\Ep) when P G {p, a} we 

must integrate equation (17) over all momenta P]r that yield Ep and make use of the fact 



that Pjf^ = ^p. 



Pp{Ep\Ep) = j d'V'p Pp{P'p\Ep) 5{Ep - E'p) (18) 

where E'p is the energy associated with P^, i.e. E'p + rriN = ^/{P^pf~+~rn^ . After a little 
effort this produces 



Pp{^Ep\Ep) 



Ep + rriN 



{Pp-Ppf\ (-{Pp + P, 



^"Pt — 24 — )-'^^\ 24 



(19) 

Likewise, when P e {Z} we can obtain Pp{Ep\Ep) but the relation between P^^ and Ep is 
now through a Lorentz transformation and consequently yields the more complicated 

Pp{Ep\Ep) = — , exp ' 



Ppy/2^ V 2a|, 
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S^exp{-xl) - exp{-xl) + ^ [erf (a;+) - erf (a;_) ]| (20) 



where 7^ = Pp/ (3pmN, clf — ^fI'^ini and x+ and X- are 

- ^^^^ 

x_ = ^pMl-PppF)-fF_ ^22) 
v2 aj? 

We have used the expressions for Pp and e|, given by Morrissey (1989), 

" '°Ap(A;,-l)' ^^"^^ 

= 8 7 Pp MeV/nucleon, (24) 

Az Pp 7p 

with eo ~ 100 MeV, Az the mass number of the nucleus to be spalled (the heavier of P and 



T) while (3p = vp/c and 7p = 1/a/1 — are the Lorentz variables of the projectile. We 
have ignored the decreases in ep that occur for values of Ep < 100 MeV/nucleon indicated 
by Stokstad (1984): at such small energies, and with our choices of A, the value of Sp is very 
close to unity and every fragment isotope is captured by the Galaxy rendering the decrease 
in ep irrelevant. The expression for Pp in equation (24) diverges as /3p becomes small so it 
must be replaced by another in this limit. The substitute we have used is 

pF^lpf3p , rriNpp, (25) 

Zip -\- Zip 

where Al is the lighter mass of the interacting nuclei P and T. This equation assumes the 
fragment carries a fraction Ap/{Ap + Ap) of the projectile momentum in the frame where 
the heavier nucleus is at rest and the changeover from equation (24) to (25) again occurs at 
energies < 100 MeV/nucleon. Gumming, Haustein & Hseuh (1981) state that such a change 
is expected to occur at roughly this energy. 

The two probability distributions are shown in figures (1) and (2) for the interaction of 
p/H and ^^O to produce ^Be when the Lorentz factor of the projectile is 7p = 2. The figures 
show that the emerging Beryllium energy is peaked below ~ 5 MeV/nucleon for the forward 
case and that the spread in Beryllium energy for the inverse reaction is considerable. 

With the probability distributions now under control it is straight forward to calculate 
the rates of production and the flux by inserting the distributions into equations (13) and 
(14). We compile the rates again for A = 5 g/cm^ and 10 g/cm^ which we show in table (2) 
again only for the reactions involving p/H with ^^O. After comparing with the results in table 
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(1) we can see that straight- ahead approximation accurately predicts the production rate 
for the forward reactions. This result was to be expected: the distribution of the fragment 
energy in localized to < 5 MeV/nucleon as shown in figure (1) and at these energies the 
range (the distance travelled by a particle before coming to rest per unit mass density of 
the medium) is much shorter than A so all fragments are trapped. After examining the 
inverse reactions we see a slight increase (< 8%) in the production rates relative to the 
straight-ahead calculations but the differences are small and do not warrant rejection of the 
approximation. The differences are of a similar magnitude as those found by Tsao et al. 
(1995). The approximation's success again has a simple explanation: the dispersion of the 
Beryllium energy ese is only of order ~ 10 MeV and so Pbc ~ Po, E^e ~ Eq. 



4. TwoStep Reactions 

So far we have only considered the first term in the reaction expansion from equation 
(7) but now we turn our attention to the second term, Q^^ ^ to determine its contribution to 
the rate of change the abundance. 

If the fragment Fi from the reaction P + Ti Fi undergoes a subsequent reaction 
F1 + T2 ^ F2 then from equation (7) we have 

POO 

Q'iliEF,) ^^^1 <f>F!iEF,) riT, T^^/'iEF,) P{Ef,\E^J (26) 

Inserting the expansion for (j)p^ and the expression for from equation (12) then 

gg(i^,J= Y: ^^Z^H rdE,,^>^-^{E,,)alf\E,,)P,,[E,,\E,^ (27) 
p,Ti,r2 Pl^H Jo 

where the flux ^^^"^^(i^pj appearing in equation (27) is the normalized onc-stcp flux that was 
deflned in equation (12). Introducing the quantity Qp^^''^^ as the contribution to equation 
(27) from each triplet P, Ti, T2 and then inserting the result into equation (10) allows us to 
define a normalized two-step Rp^^'^'^ by 

poo 

i^^'^^ = E / dEF. QW"^'iEF.) SfAEf.) (28) 
P/nH ^ Jo 

so that 

J (2) 
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We have calculated these two-step reaction rates R^^ ^' ^ again for A = 5 g/cm^ and 10 g/cm^ 
using our improvement of the fragment energy distribution from section §3 and in table (3) 
we show the results for the reactions involving p/H and ^^O. For the forward cases, i.e. 
when P G {p,a}, we find that the two-step rates vary in magnitude from ~ 10^^^ mgcm^ 
for ^Li to ~ 10"^'' mgcm^ for ^^B whereas the inverse rates, i.e. P G {Z}, are all around 
~ 10^"^^ mgcm^. The forward rates are very sensitive to the threshold energies of the 
reactions F1 + T2 — > F2 since the fluxes of Fi peak at very low energies as indicated by figure 

(1) . In both cases these rates certainly appear much smaller than those in tables (1) and 

(2) and it is tempting to regard the two-step rates are negligible. However this would be 
an erroneous conclusion because attention must be given to the pre-factors of equation (28). 
For an ISM of only H and He with a mass fraction of He of Y the p/nn term is 



Including this factor we see that the forward rates arc at least five orders of magnitude 
smaller than the one-step rates in table (2). In contrast, when we apply the p/nn prc-factor 
to the inverse rates that we find they arc only ~ 1/10 smaller! This same result was found 
by Ramaty et al. (1997) who also saw substantial changes to their calculations of the ratios 
of production rates when they included the two-step contributions. 

This result may seem surprising but it is not entirely unexpected. The inefficient retar- 
dation of the CRs means that the flux of intermediaries Fi accumulates over a long period 
of time. The flux of Fi is only ~ 10^^ — 10^'^ smaller than the primary (jy'p . After we 
factor in the smaller average energy of the intermediaries (which will result in a higher frac- 
tion of trapped fragments) and the multiplicity of the two-step channels we obtain two-step 
production rates that are comparable with the one-step. 

The signiflcant two-step contribution to the rate of change of abundances means that 
we must evaluate the higher order terms in the expansion of Q^p . At the two-step level the 
number of reaction rates we must compute is considerable, if we find equally large contri- 
butions from the three-step reactions and greater then the computational burden becomes 
excessive. At this point we would have to abandon this approach to spallation calculations, 
so we must either show that the three-step contribution is negligible or we must find an 
alternative and more efficient approach to the calculation of production rates. 



1 



l-Y 



(30) 



p/riH 



rriH 



5. 



A Cascade Calculation 



The three-step rates are, in principle, simple to compute since we proceed in the same 
fashion as the one and two-step equations we derived in sections §2 and §4. In practice 
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there are so many pathways between any given primary projectile Pq and the fragment F3 
under consideration that the procedure becomes increasingly laborious. Instead we adopt the 
approach of Moskalenko et al. (2002) and utilize the simple property that in any spallation 
reaction the fragment is always lighter than the heavier reactant. 

We begin by splitting the flux 0^ of any nucleus A into the sum of two components, the 
primary accelerated from the ISM and then the secondary flux The secondary flux 
0^'' is dependent upon the fluxes of all the nuclei B heavier than A which again are sum of 
primary and secondary components. Thus we have (f)A = 4>a^ + 4>a^ and 0^'' = X1b>a /(*i^B'' + 
0^'') where / represents the successive application of the linear equations (5) and (1). The 
secondary flux of the lightest member of B is expressible as 0^'' = "^cyB fi'P^c + '^c'') ^"^^ 
so it can be eliminated from the function for (j)'^ . We continue consecutive elimination of the 
lightest secondary flux in the expression for 0^-* until we reach the heaviest nucleus we are 
considering. The flux of this heaviest nucleus does not have a secondary flux by construction 
so we flnd that is a function only of primary fluxes 0^°^. We therefore introduce the 
quantity 0^° as the contribution to from each primary flux which we label by Pq 



where we use the symbol F' instead of P to emphasize that the second term involves reactions 
between fragments and the ISM. The set of fragments F' is a subset of the Z nuclei so the 
ISM targets T in the second term of equation (32) can only be H or He (we have been ignoring 
Z -\- Z' reactions) and therefore the fragment F must be lighter than F' . We introduce the 
quantities as the contribution to this equation from each term of the expansion over Pq 
and 'normalize' the primary flux as ^^j* = z/'pg in exactly the same way as section §2 in 
order to extract the dependence upon the composition of the CRs. After their introduction 
we reach the most important expression of this section: 




(31) 



Inserting this expression into equation (5) for 0p we obtain 




(32) 




+ J2 J2 dEF>(t)l%EF>)nTap'^{EF.)PF{EF\EF,). 




(33) 
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We have introduced the symbol F + 1 as the hghtest fragment F' from which F may be 
produced in the reaction with T while the upper limit to the sum, Pq — 1, is the heaviest 
nucleus that can be formed after an interaction of Pq with any component of the ISM. We 
have denoted this nucleus suggestively by Pq ~ 1 but we note that in the case of Pq G {p, «} 
the heaviest intermediary F' is formed from the heaviest component of the ISM. Equation 
(33) shows that the source Qp° is the sum of the one-step contributions Pq + T — > P plus 
the reactions involving all secondary nuclei F' and T. This equation apphes to any fragment 
F and as the mass of F increases the number of intermediaries becomes smaller. However a 
major change occurs when F = Pq — 1 because for this nucleus there are no intermediaries 
that can produce P, only the one-step process involving Pq and T manufacture F = Pq — 1. 
For this nucleus the source spectrum Q^=p^_i is calculated with only the first half of equation 
(33) i.e 

poo 

Y^VTUH / dPpo^Po(PpJa;°^(PpjPp„_i(Pp,_i|PpJ. (34) 

rp Jo 

The spectrum Q^Lpo-u hence the flux (f)p°^p^^__-^^, are proportional to ap,j and therefore 
we introduce the normalized flux in the same fashion as we did in section §2 so that 

<t>F.p,-i{EF.p,-i) = ap.*Ftf..-i(£^F-P.-i) 

Wpo-l[r^Po-l) P JEp^_i ^Po-i[^Po-i) 

Every quantity on the right-hand side of equation (35) is known so the flux ^pLpo-i is exact. 
Now that we have 0p°_i we turn our attention to QpLpo-2' heaviest fragment that can 
be produced from P' = Pq — 1. For F — Pq — 2 the sum over intermediaries P' in equation 
(33) involves only the fragment P' = Pq — 1 and we insert our solution for from 
equation (35). Both and (t>^i=z-i already known so the source spectrum (5pLz-2' 
and hence the flux 0p°=2_2, may also be calculated exactly. As we step down through the 
nuclei the source spectrum Q^" of any fragment P is only dependent upon fluxes that are 
previously derived and so we do not introduce any errors into Even the a — a fusion 

reactions that produce ^Li and ^Li can be accommodated in this scheme by first calculating 
and then appropriately inserting this solution into the expression for QeLi- For every 
fragment P the source spectrum is proportional to ap^ and so we write the rate of change 
of abundance for P for a given primary as 

^ = ap, = dE'pQ^°{E'p) Sp(E'p) (36) 

which defines the normahzed cascade rate Rp" . 
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Previously we extracted the dependence upon the ISM abundances j/t but we have 
refrained in this cascade calculation because if we do then we run into the same problem as the 
reaction expansion, namely, vast arrays of rates to calculate for all the possible pathways from 
the primary Pq to F. In order to calculate the source spectra, the fluxes of the intermediaries 
and the rate Rp' we need to specify a composition for the ISM. This introduces a chemical 
dependence we have been at pains to avoid but the sufficient accuracy of the one-step rates for 
the forward reactions, when Pq G {p,a}, means that we only follow this cascade for Pq e {Z} 
which is why we called the heaviest intermediary Pc, — 1. In the inverse calculations the heavy- 
element content of the ISM is irrelevant since all the reactions are with H and He as targets 
so we only need specify the abundance of Hehum in the ISM. The dependence upon the 
ISM chemistry is therefore very weak since yue is essentially constant. In table (4) we show 
the cascade rates Rp' for the case when the primary projectile is Oxygen together with the 
rate as calculated by using the reaction expansion up to the two-step order. The grammage 
is set at the constant value of A = 5 g/cm^, we use the straight- ahead approximation and 
the helium abundance is set to yuc = 0.1. It is clear from the table that the cascade and 
reaction-expansion rates agree to < 3% and thus we conclude that terminating the reaction 
expansion at the two-step level gives sufficiently accurate rates that the difference between 
the two methods is neghgible. 



6. Conclusions 



The results of any calculation of the spallogenic rates of production for Lithium, Beryl- 
lium and Boron often find their application in trying to remove this synthesis mechanism 
from the observations in order to elucidate the other contributions to the abundances. Sim- 
ilarly, the ratio of the production rates has been frequently used to infer the spallogenic 
abundance of one of the LiBeB elements after observation of another. Both applications rely 
on the accuracy of the calculated rates so erroneous results can occur if the estimates for the 
contribution to the abundance from spallation are wrong. There are many aspects of the 
calculation where errors can enter but two do not have astrophysical origins: the straight- 
ahead approximation and the reaction expansion. We find that relaxing the straight-ahead 
approximation does increase the complexity of the problem because the distribution of the 
fragment energies must be taken into account but that using a more sophisticated description 
produces changes that are less than ~ 5%. This is too small to insist that the straight-ahead 
approximation be rejected at the present time. Much larger changes were found when we in- 
cluded the two-step terms from the reaction expansion so that our faith in this methodology 
waned. However, after computing the rates to all order we found changes that the difference 
between the cascade and two-step rates is only a few percent and therefore sufficiently ac- 
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curate rates can be calculated by terminating the reaction expansion at the two-step term. 
We therefore conclude that the errors arising from two assumptions in the nuclear physics 
aspects of the calculation can be removed leaving only that from the uncertainty in the 
cross-sections. 

The authors would like to thank Gary Steigman for useful discussions. We acknowledge 
the support of the DOE through grant DE-FG02-91ER40690. 
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Table 1. The production rates Rp using the straight-ahead approximation for the case 
when the reactants are p/H and ^^O. We have used two, different, constant values for the 
grammage of A = 5 g/cm^ or A = 10 g/cm^. The reaction rates are in units of cm^. 



P, T, F A = 5 g/cm^ A = 10 g/cm^ 



P, ^^0, 


6Li 


1.29 


X 10" 


-26 


1.29 


X 10" 


-26 


p, 1^0, 


^Li 


2.06 


X 10- 


-26 


2.06 


X 10" 


-26 


P, ^^0, 


9Be 


4.12 


X 10- 


-27 


4.12 


X 10- 


-27 


p, 1^0, 




1.49 


X 10- 


-26 


1.49 


X 10" 


-26 


P, ^^0, 




2.80 


X 10- 


-26 


2.80 


X 10" 


-26 


1^0, H, 


6Li 


2.03 


X 10- 


-27 


2.94 


X 10" 


-27 


1^0, H, 


^Li 


3.45 


X 10" 


-27 


4.97 


X 10- 


-27 


i«0, H, 




6.77 


X 10" 


-28 


9.87 


X 10" 


-28 


i«0, H, 




3.14 


X 10" 


-27 


4.37 


X 10" 


-27 


i«0. H. 




7.55 


X IQ- 


-27 


9.86 


X IQ- 


-27 
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Table 2. The production rates Rp calculated after relaxing the straight- ahead 
approximation for the case when the reactants are p/H and ^^O. As in table (1) we present 
the results for two, different, constant values for the grammagc: A = 5 g/cm^ and 
A = 10 g/cm^. The reaction rates are in units of cm^. 



P,T,F 


A = 


5 g/cm^ 


A = 10 g/cm^ 


p, 1^0, %i 


1.29 


X 10^ 


-26 


1.29 


X 10" 


-26 


p, 1^0, ^Li 


2.06 


X 10- 


-26 


2.06 


X 10- 


-26 


p, 1^0, ^Be 


4.12 


X 10" 


-27 


4.12 


X 10- 


-27 


p, 1^0, lOB 


1.49 


X 10- 


-26 


1.49 


X 10- 


-26 


p, 1^0, iiB 


2.80 


X 10" 


-26 


2.80 


X 10- 


-26 


1^0, H, "Li 


2.19 


X 10- 


-27 


3.09 


X 10- 


-27 


1^0, H, ^Li 


3.67 


X 10- 


-27 


5.16 


X 10- 


-27 


1^0, H, ^Be 


7.17 


X 10- 


-28 


1.02 


X 10- 


-28 


1^0, H, lOB 


3.24 


X 10" 


-27 


4.45 


X 10- 


-27 


1^0, H, i^B 


7.73 


X 10- 


-27 


9.99 


X 10- 


-27 
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Table 3. The production rates Rp ^' ^ calculated without using the straight- ahead 
approximation for the case when all the reactants are p/H and ^^O. As in table (1) we 
present the results for two, different, constant values for the grammage: A = 5 g/cm^ and 
A = 10 g/cm^. The reaction rates are in units of mgcm^. 



P, Ti, Ts, F A = 5 g/cm^ A = 10 g/cm^ 



p, 1^0, H, 


6Li 


5.98 


X 10" 


-53 


5.99 


X 10- 


-53 


p, 160, H, 


^Li 


2.44 


X 10- 


-53 


2.44 


X 10" 


-53 


p, 160, H, 


"Be 


1.74 


X 10" 


-57 


1.74 


X 10- 


-57 


p, 160, H, 


log 


2.23 


X 10- 


-61 


2.23 


X 10- 


-61 


p, 160, H, 


iiB 


1.50 


X 10" 


-67 


1.50 


X 10- 


-67 


160, H, H, 


6Li 


1.04 


X 10- 


-48 


2.60 


X 10- 


-48 


160, H, H, 


^Li 


1.23 


X 10- 


-48 


3.01 


X 10- 


-48 


160, H, H, 


9Bc 


7.29 


X 10- 


-49 


1.80 


X 10- 


-48 


160, H, H, 




1.76 


X 10- 


-48 


4.27 


X 10- 


-48 


160, H, H, 


iiB 


2.92 


X 10" 


-48 


7.14 


X 10- 


-48 
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Table 4. The production rates Rp° for the case when the primary is ^^O and the rate as 
determined by using the reaction expansion up to the two-step terms. The path length was 
set at the constant value of A = 5 g/cm^ and we used the straight-ahead approximation. 

The reaction rates are in units of cm^. 



^0, 


F 


cascade 




reaction expansion 


1^0, 


6Li 


2.79 X 10- 


-27 


2.72 X 10-27 


1^0, 


^Li 


4.38 X 10" 


-27 


4.33 X 10-27 


1^0, 


9Be 


1.19 X 10- 


-27 


1.15 X 10-27 


1^0, 




4.50 X 10" 


-27 


4.42 X 10-27 


160, 




9.79 X 10- 


-27 


9.77 X 10-27 
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Fig. 1. — The probability distribution P{Ef\Ep) as a function of the fragment energy Ep 
for the reaction H O — >^ Be at 7jj = 2 and A = 5 g/cm^. 
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0.012 



0.01 - 




1.5 



1.75 



2 
Yo 



2.25 



2.5 



Fig. 2. — The probability distribution P{Ef\Ep) as a function of the Lorentz factor 'fp of 
the fragment for the reaction ^^O + H — >^ Be at 70 = 2 and A = 5 g/cm^. 



